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TECHNICAL MEMORANDUM X-64662 

CONJUGATE GRADIENT DETERMINATION OF OPTIMAL 
PLANE CHANGES FOR A CLASS OF THREE-IMPULSE 
TRANSFERS BETWEEN NONCOPLANAR CIRCULAR ORBITS 

INTRODUCTION 

Recently the problem of optimally distributing three plane changes 
when impulsively transferring between two noneoplanar, circular orbits was 
reconsidered. The particular version studied involves a velocity impulse at 
the inner circular orbit that places the vehicle at the perigee of an ellipse 
whose apogee is greater than the radius of the outer circular orbit. At the 
apogee of this ellipse, another velocity impulse places the vehicle on an 
ellipse whose perigee coincides with the outer circular orbit. At the perigee 
of this ellipse, another velocity impulse occurs to place the vehicle in the 
outer circular orbit. Each of the three impulses can involve a plane change, 
and the problem is to minimize the total velocity impulse, i.e. , the sum of 
the velocity impulses, by defining the optimal plane change at each impulse. 
This problem is interesting from two standpoints. First, when the radii of 
the two circular orbits and the transfer ellipse apogee are close in magnitude, 
obtaining numerical solutions with current computer programs which are 
demonstrably very good programs was very difficult, if not impossible, for 
reasons discussed later. Secondly, the near-earth on-orbit maneuvering 
done by the Space Shuttle vehicle or Space Tug could involve just exactly the 
conditions leading to numerical difficulties; therefore, it is important for 
planning and other purposes that accurate numerical results be obtainable 
and available. The first point primarily instigated this report, and led to 
the development of an extremely good numerical algorithm. 


PROBLEM DESCRIPTION 


To fix ideas, consider the sketch in Figure 1 of the geometry involved. 
At point 1, sufficient velocity is imparted to the vehicle to achieve perigee 
velocity of the first half of the transfer ellipse. At point 2, velocity is added 
to obtain apogee velocity of the second half of the transfer ellipse. At point 3, 
the vehicle velocity is circularized. The following formulas define the various 
velocities: 
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Figure 1. Geometry of 
the orbital transfer. 
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where the subscripts c, a, and p refer 
to circular, apogee, and perigee, 
respectively, and m is the gravitational 
parameter of the earth. 


The velocity impulses in question are 
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where V - and v a + refer to the apogee velocities of the first half and 
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second half of the transfer ellipse. Including the possibility of noncollinear 
velocity impulses, the AV* s may be written as 
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where ^ defines the angle between the plane of the initial, inner circular 

orbit and the plane of the first half of the transfer ellipse; & defines the angle 

between the planes of the first half of the transfer ellipse and the second half 
of the transfer ellipse; and ^ defines the angle between the plane of the second 
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half of the transfer ellipse and the plane of the outer circular orbit. Factoring 

out — and performing other manipulations allow these equations to be written: 
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The total Ay is 


AV total ’ AV 1 + AV 2 + AV 3 • 
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and if it is referenced to the initial circular velocity, /— , AV, , . may be 
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Making the following replacements: 
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the nondimensional (normalized) total impulse (also referred to as the 
payoff) is 
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AV n 


2 V 2 1/ 

= (i + H 1 - 2H 1 cos 1 ? 1 ) /2 + H 2 (i + H 3 - 2H. 3 cos,? 2 ) /z 
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The problem statement reads as follows: Choose £ , i? , and to 

f 1 m u 

minimize Av,^ subject to the equality constraint 

*1 + ^2 + ^3 = ’ (2) 


where ii^is the angle between the planes of the initial and final circular orbits, 
and the inequality constraints are 

t>i s 0, i> 2 2: 0, ^ > 0 • (3) 


Reasons based on physical grounds indicate that the inequality constraints 
can be neglected in the mathematical solution. However, in a later section, 
it is demonstrated how even more restrictive inequality constraints are 
treated and solved. The equality constraint is easily accounted for by solving 
equation (2) for £ , for example, and substituting into equation (1). The 
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minimization then proceeds with only two free variables, and ■& . Since 

t 12 

the free variables enter into AV in the arguments of trigonometric functions. 


a numerical procedure is indicated for the minimization process. Invariably 
the gradient of the payoff will enter into the minimization process either to 
establish a descent direction or to indicate when a minimum has been achieved, 

since a necessary condition for a minimum is that the gradient equal zero. 
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The gradient of AV^ is 
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Comparing equations (i) and (4) shows that the same radicals are involved 
if; ) in each. Previous experience showed that when the apogee of the transfer 
ellipse was much larger than the radius of the inner circular orbit, the 
numerical minimization proceeds without incident. : However,, for the close - 
in'orblFtransferS^ ^ considered herein, one A an^ possibly, two. d’ifficult'ies . apj.se: 

1 ) ’ the' first evolves because the radicalsofequations (1) .and.. (4) . involve 
the aiffe renc ing of quantitie s , each of which is about two, . in magnitude , apd 
mis subtraction results in a serious loss of significant digits; and (2), if the 

^ ■ . - - ••*/* -'-I . • / i'j?, y- ^ r tri TiV : > "Vli „ V ' ' ■ * •• r ‘ i : 

numerical minimization procedure is a zero finding algorithm wherein the 
^^minimum is^ defined by finding the roots of equation. (4) ; the Jacobian of ~i 
equation (4) can be involved. As it re suits , for the close-in transfers con- 
- Brio • red here , two' of the angles are very smaHi ‘In terms of the small 

y Xddifiference. between the cosines of these small angles, and one, : it can.be shown 
/that; terms appear :in the Jacobian whose order.of magnitude is .the square-i 
b*v, ^qf these differences. Typically, for differences of the. order 10 


-3 


, terms ob- 
This extreme 


tained by -subtraction appear in the Jacobian- of the order of 10 
loss of precision dictates that most standard gradient or Newton-Raphson 
techniques will behave erratically and, when they do converge, will converge 
to false minima. This type of performance was in fact observed with some 
very good existing computer programs. 
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The answers to both difficulties lie in a detailed inquiry into and analy- 
sis of the underlying conditions causing the loss of precision. This task is 
simplified if an analysis of the Jacobian can be eliminated and consideration can 
be limited to the payoff and its gradient. This is the case for conjugate gradi- 
ent numerical minimization techniques of which a Sorensen [1 ] modified 
Fletcher-Reeves [2 ] version was chosen for the numerical results reported 
here. The version used here will be discussed and compared extensively 
with other techniques in a separate report. Let it suffice to say that conjugate 
gradient algorithms proceed along successive paths of descent with local 
minima of the payoff being found along each path. At each iteration, only the 
payoff and its gradient need be known. Assuming accurate numerical infor- 
mation, successive iterations are guaranteed to decrease the payoff until 
final convergence is achieved. A basic property is that a payoff quadratic 
in n variables will be minimized in n iterations. 

NUMERICAL ANALYSIS 

The numerical analysis begins by observing that H,, H , and H are 

J. o o 

all greater than one for the problem considered. Taking as typical, let 
H A = 1 + A 1 . (5) 

Now, 
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Substituting equation (5) into the right side of the preceding equation yields 

2 V 2 2 ”1 V 

( i + - 2H^ cos 72 = (i + A^ - cos i^) + sin 72 . 


Further, let cos ^ = i - 6^; then, 
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Equations (6) are accurately computable if A. and 6. are accurately available. 
To obtain A., consider as typical again. From its definition, 
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Since a < 1, the binomial series converges and represents H^: 
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r - r 
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All these series are valid for ia.l < 1. 


The Maclaurin series for cos £ is 


2 4 

cos = 1 - — + - 4f + ••• ("D 
-t* • • • n — 1,2,3,.,. , 


2(n-l) 
(2n-2) 1 


so that 


■>, 2 ,> 2n 
6, - i - - -jj- - — + .. (-I) n ^ 


• • • n “ lj 2, 3, • • • j 


(8a) 


6 3 = 


_ 2 _ 

2 ! 


2 

41 


+ ... (- 1 } 


n 2 


2n 


(2n) ! 


. n = 1, 2, 3, . . . 


and 


(8b) 


i? 


6 5 - 


3_ 

2! 


3 

4! 


2n 


+ . 


. (-D n 3 


(2n) 1 ® 1> 2, 3, . . » 


(8c) 


These series are valid for all finite 

The series in equations (7) and (8) are alternating series so that the 
error committed by truncating them at any term is less than the value of the 
first term in the remainder. Assuming that equation (7) will be used only 
if | I < 0.1, an estimate of the number of terms required to limit the 
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error to < 10“ 18 is obtained as follows. It is easily seen that the magnitude 

of each coefficient, a , inequation (7) is =£ — . Thus, 

n 2 n 

la | (0.i) n < -°‘ — } - = i0 _i8 
n 2 n 

or 

-n - n log 2 =• -18; 
however, 

log 2 fs 0.3 
Therefore, 

1.3n = 18 
or 

n « 14. 

Therefore, 13 terms of the series in equation (7) are more than sufficient 

”18 

to compute A. with an error less than 10 when a. < .1. The magnitude 

u i i 2n 

of the coefficients inequation (8) is , . Limiting |t? . | < 2 

(w 114. 6°), an error estimate for 13 terms is 


2 28 

(28)1 


2.7 x 10 
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3 x 10 


29 


= 0.9 x 10 


-21 


~18 

That is, 13 terms are more than sufficient to limit the error to less than 10 
No effort was made to sharpen these error estimates or economize the series 
used. 
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UNCONSTRAINED NUMERICAL EXAMPLES 


The foregoing simple analysis is sufficient to eliminate the difficulties 
discussed earlier. Interestingly, numerical results indicate multiple solutions, 
the solution obtained being dependent on the starting point since a gradient 
technique will, or at least should, go to the bottom of whatever valley it begins 
in. Recalling the implicit inequality constraints of equation (3) , the solution 
points lie on a plane as shown in Figure 2. The payoff along the boundary of 
this plane can be explicitly computed. As it results the true minima lie on 

the interior of this plane. These were 
found by beginning the algorithm at 
each of the vertices in turn and proceed- 
ing to the local minimum. All of the 
minima lie close to a vertex, and only 
a maximum of three was ever found, 
even though in some cases many starting 
points within the plane were attempted. 
The global minimum is found by compar- 
ison of the local minima. It is not the 
purpose of this report to present an 
analysis of close-in orbital transfer; 
therefore, Tables 1 through 8 are pre- 
sented only to show the speed and accu- 
racy with which solutions are obtained 
for representative situations. A general 
feature of these tables is that the global 
minimum occurs when the maximum plane change takes place at the maximum 
radius, i.e., apogee, of the transfer ellipse. 



CONSTRAINED NUMERICAL EXAMPLES 


An easily applied transformation technique can be used to account 
for even more general inequality constraints than those in equation (3) . 
Consider inequalities of the form 

L. < < U. i = 1,2,3 

ill 

where L. and U . are lower and upper bounds respectively on & .. Then, the 
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following transformations automatically satisfy these inequalities : 

, X 2 

i?. = L. + (U. - LJ sin X- , 

11 11 i 

where y. becomes the new minimization variables. These transformations 
i 

need to be carefully applied. To avoid an unnecessarily complicated expla- 
nation, the application used here is to consider explicit limitations on and 
i ?2 °f the form 


= c sin x 1 


and 
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where c > 0 and c > 0 and e + c ^ # 
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These introduce the partial s - — • and - — 

sx t sx 2 


into the gradient equations (4) 


in a simple way; i.e. , 
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dX A 
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dX c 


The constrained minimization results depend on the values given c . and c , 

and it might be expected that either, both, or neither t? and S- lie on their 

1 2 

respective boundaries. Table 9, which can be compared to Table i, illu- 
strates all four possibilities. It shows that this type of constrained minimi- 
zation is no more difficult than unconstrained minimization. 


CONCLUSIONS 


The numerical difficulties of determining optimal plane changes associ- 
ated with a particular class of three -impulse transfer between noncoplanar. 
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circular orbits have been eliminated using simple series. The conjugate 
gradient algorithm developed resulted in an ideal tool to obtain very accurate 
solutions, simply and rapidly, for unconstrained and constrained angles. 


TABLE 1. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 150 N. MI. CIRCULAR 
ORBITS INCLINED AT 28.5 DEGREES (TRANSFER APOGEE IS 200 N. MI.) 


Iteration 


*2 

3 

Gradient 

Payoff 

0 

0.0 

0.0 

28.5 

0.56E-03 

0.50013379 

1 

1.46114 

1.46114 

25.5777 

0. 15E-08 

0.49216636 

2 

1.51796 

1.40452 

25.5775 

0.43E-09 

0.49216479 

3 

1.49644 

1.33453 

25.6690 

0.53E-11 

0.49216410 

4 

1.49324 

1.33551 

25.6712 

0. HE -12 

0.49216410 

5 

1.49344 

1.33683 

25.6697 

0.50E-18 

0.49216410 

6 

1.49344 

1.33683 

25.6697 

0.0 

0.49216410 

0 

0.0 

28.5 

0.0 

0.28E-03 

0.49333864 

1 

0.0 

27. 8431 

0.656886 

0.28E-03 

0.49091133 

2 

1.22383 

26.6175 

0.658645 

0.20E-08 

0.48613730 

3 

1.22387 

26.6016 

0.674486 

0.20E-12 

0.48613706 

4 

1.22401 

26.6016 

0.674347 

0.73E-13 

0.48613706 

5 

1.22 423 

26.6013 

0.674492 

0.72E-19 

0.48613706 


28. 5 

0.0 

0.0 

0.29E-03 

0.50096085 

1 

27. 5403 

0.0 

0.959736 

0.29E-03 

0.49823819 

2 

1.23454 

26.2510 

1.01448 

0.27E-06 

0.48621170 

3 

1.43016 

26.2553 

0. 814521 

0.39E-07 

0.48617471 

4 

1.42698 

26.3974 

0.675637 

0.44E-07 

0.48615935 

5 

1.22422 

26.6023 

0.673492 

0.74E-11 

0.48613706 

6 

1.22423 

26.6013 

0.674492 

0.11E-16 

0. 48613706 

7 

1.22423 

26.6013 

0.674493 

0.39E-17 

0.48613706 

8 

1.22423 

26.6013 

0.674492 

0.43E-23 

0.48613706 


a. Gradient 
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14 















TABLE 2. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 150 N. MI. CIRCULAR 
ORBITS INCLINED AT 60 DEGREES (TRANSFER APOGEE IS 200 N. MlJ 


Iteration 


*2 

*3 

a 

Gradient 

Payoff 

0 

0.0 

0.0 

60.0 

0.45E-03 

1.0051436 

1 

0.546747 

0.546747 

58.9065; 

0.17E-05 

1.0000558 

2 

0.707780 

0.404563 

58.8877 

0.61E-07 

0.99991549 

3 

0.688537 

0.376323 

58. 9351 

0.63E-08 

0.99991119 

4 

0.678625 

0.383077 

58.9383 

0.34E-09 

0.99991071 

5 

0.676701 

0.381311 

58.9420 

0.10E-12 

0.99991069 

6 

0.676738 

0.381284 

58.9420 

0.14E-18 

0.99991069 

7 

0.676738 

0.381284 

58.9420 

0.50E-20 

0.99991069 

0 

0.0 

60.0 

0.0 

0.22E-03 

0.99138951 

1 

0.0 

59.6739 

0.326106 

0.22E-03 

0. 98974128 

2 ' 

0.641237 

59.0316 

0.327157 

0.24E-08 

0. 98646520 

3 

0.641244 

59.0287 

0.330028 

0.24E-14 

0. 98646515 

4 

0.641246 

59.0287 

0.330026 

0. 10E-14 

0. 98646515 

5 

0.641251 

59. 0287 

0.330028 

0.44E-21 

0. 96640515 

0 

60.0 

0.0 

0.0 

0.23E-03 

1.0103875 

1 

59.6427 

0.0 

0.357284 

0.23E-03 

1.0086236 

2 

59.2477 

0.394300 

0.357992 

0.84E-09 

1.0067450 

3 

59.2456 

0.394303 

0.360086 

0.68E-14 

1.0067449 

4 

59.2456 

0.394308 

0.360082 

0.20E-14 

1.0067449 

5 

59.2456 

0.394313 

0.360086 

0.44E-21 

1.0067449 


a. 


Gradient = 



+ 
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TABLE 3. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 150 N. MI. CIRCULAR 
ORBITS INCLINED AT 28. 5 DEGREES (TRANSFER APOGEE IS 151 N. MI.) 


Iteration 

J 

*1 

^2 

^3 

. Gradient 3- 

Payoff, . 

— 

0 

0.0 

0.0 

28.5 

0. 56E-03 

0.49593920 

1 

0.799631 

0. 799631 

26.9007 

0.14E-07 

0.49066539 

2 

0. 734489 

0. 865430 

26. 9001 

0.23E-09 

0.49065984 

3 

0. 742474 

0.875737 

26.8818 

0.36E-11 

0.49065974 

4 

0.741365 

0. 876595 

26.8820 

0.17E-12 

0.49065974 

5 

0.741536 

0.876955 

26.8815 

0.55E-18 

0.49065974 

6 

0.741536 

0.876956 

26. 8815 

0.40E-26 

0.49065974 

0 

0.0 

28.5 

■ 0.0 

0.28E-03 

0.49155631 

1 

0.0 

28.4847 

0.015318 

0.28E-03 

0.49150473 

2 

0.703079 

27.7816 

0.0153418 

0.54E-09 

0.48894035 

3 

0. 703080 

27.7813 

0.0155842 

0. 18E-16 

0.48894035 

4 

0.703080 

27.7813 

0.0155842 

0.17E-16 

0.48894035 

5 

0.703082 

27.7813 

0.0155842 

0.29E-21 

0.48894035 

0 

28. 5 

0.0 

0.0 

0.29E-03 

0.49675259 

1 

28.4811 

0.0 

0.0189021 

0.29E-03 

0.49669797 

2 

27.4330 

1.04801 

0.0189459 

0.12E-08 

0.49389062 

3 

27.4324 

1.04801 

0.0196373 

0.50E-15 

0.49389061 

4 

27. 432 4 

1.04801 

0.0196367 

0.48E-15 

0.49389061 

5 

27.4323 

1. 04805 

0.0196374 

0. 50E-23 

0.49389061 





a. Gradient 








TABLE 4. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 150 N. MI. CIRCULAR 
ORBITS INCLINED AT 60 DEGREES (TRANSFER APOGEE IS 151 N. MI. ) 


Iteration 

*1 

*2 

^3 

Gradient 3, 

Payoff 

0 

0.0 

0. 0 

60.0 

0. 45E -03 

1.0001036 

1 

0.351493 

0.351493 

59.2970 . 

0.66E-08 

0-99659350 

2 

0.346132 

0.356896 

59.2970 

0.47E-11 

0.99659319 

3 

0.3462 73 

0.357045 

59.2967 

0.14E-14 

0.99659319 

4 

0.346270 

0.357047 

59.2967 

0.31E-17 

0.99659319 

5 

0.346270 

0.357047 

59.2967 

0.73E-25 

0.99659319 


0.0 

60.0 

0.0 

0.22E-03 

0. 99470577 

1 

0.0 

59.9931 

0.00687260 

0.22E-03 

0.99467142 

2 

0.341285 

59.6518 

0.00688441 

0.14E-09 

0. 99293887 

3 

0.341285 

59.6518 

0.00689966 

0.88E-19 

0.99293887 

0 

60.0 

0.0 

0.0 

0.23E-03 

1.0053444 

1 

59.9928 

0.0 

0.00716288 

0.23E-03 

1.0053089 

2 

59.6249 

0.367966 

0.00717615 

0.83E-09 

1.0035026 

3 

59. 6248 

0.367966 

0.00721726 

0. HE -18 

1.0035026 

4 

59.6248 

0.367966 

0.00721726 

0.12E-18 

1.0035026 

5 

59.6248 

0.367966 

0.00721726 

0.30E-19 

1.0035026 



a. Gradient 
















TABLE 5. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 110 N. MI. CIRCULAR 
ORBITS INCLINED AT 28.5 DEGREES (TRANSFER APOGEE IS 150 N. MI.) 



*1 

2 

's' 

Gradient 3, 

Payoff 

0.0 

0.0 

28.5 

0.57E-03 

0.49650885 

0.638691 

0.638691 

27.2226 

0. 15E-06 

0.49342843 

0.975791 

0.312526 

27.2117 

0.55E-07 

0. 49334451 

0. 947248 

0.197258 

27.3555 

0.23E-07 

0. 49332974 

0.807265 

0.231993 

27. 4607 

0.43E-08 

0.49331817 

0.789268 

0.208230 

27. 5025 

0.46E-10 

0.49331723 

0.793417 

0.206754 

27. 4998 

0.24E-12 

0.49331722 

0.793476 

0.206919 

27.4996 

0. HE -13 

0.49331722 

0.793554 

0.206908 

27.4995 

0.34E-20 

0.49331722 

0.0 

28.5 

0.0 

0.28E-03 

0.49361548 

0.0 

27.9421 

0.557934 

0.28E-03 

0. 49160491 

0.689414 

27.2518 

0.558777 

0.27E-09 

0.48909007 

0.689421 

27.2465 

0.564077 

0.25E-13 

0.48909004 

0.689460 

27.2465 

0.564038 

0.72E-14 

0.48909004 

0.689507 

27.2464 

0.564078 

0.18E-21 

0.48909004 

28.5 

0.0 

0.0 

0.29E-03 

0.496674 

27.8294 

0.0 

0.670588 

0.29E-03 

0.49454219 

27.6133 

0.215837 

0.670905 

0.42E-11 

0.49397078 

27.6122 

0.215837 

0.672001 

0. 13E-14 

0.49397078 

27.6122 

0.215846 

0.671992 

0.15E-15 

0.49397078 

27.6121 

0.215851 

0.672001 

0.13E-21 

0.49397078 


a. Gradient = 


OAVrj 
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TABLE 6. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 110 N. MI. CIRCULAR 
ORBITS INCLINED AT 60 DEGREES (TRANSFER APOGEE IS 150 N. MI.) 



3 0.399695 0.0697055. 

4 0.367167 0.0848975 

5 0.344792 0.0753710 


6 0.346985 0.0740348 

7 0.347022 0.0740964 

8 0.347115 0.0740836 

9 0.347115 0.0740836 


0 0.0 60.0 

1 0.0 59.7323 


2 0.334814 59.3970 

3 0.334815 59.3957 

4 0.334818 59.3957 

5 0.334820 59.3957 


0 60.0 0.0 

1 59.7194 0.0 

2 59.6447 0.0745448 

3 59. 6444 0. 0745449 

4 59.6444 0.0745451 

5 59.6444 0.0745452 


1 

^3 

r -1 

Gradient 

Payoff 

60.0 

0.46E-03 

1.0041913 

59.5327 

0. 70E-05 

1.0023947 

59.4935 

0. 86E-06 

1.0021102 

59.5306 

0.28E-06 

1.0020961 

59.5479 

0.23E-06 

1.0020868 

59.5798 

0.42E-08 

1.0020823 

59.5790 

0.72E-11 

1.0020822 

59.5789 

0.14E-11 

1.0020822 

59.5788 

0.13E-17 

1.0020822 

59.5788 

0.81E-26 

1. 0020822 

0.0 

0.22E-03 

0. 99617013 

0.267709 

0.22E-03 

0.99481282 

0.268159 

0.73E-09 

0.99311326 

0.269472 

0.28E-14 

0. 99311324 

0.269470 

0.86E-15 

0.99311324 

0.269472 

0. 86E-23 

0.99311324 

0.0 

0.23E-03 

1.0052485 

0.280566 

0.23E-03 

1.0038420 

0.280749 

0.43E-10 

1.0034762 

0.281102 

0.21E-15 

1.0034762 

0.281102 

0.10E-16 

1.0034762 

0.281102 

0.20E-22 

1.0034762 
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TABLE 7. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 110 N. MI. CIRCULAR 
ORBITS INCLINED AT 28.5 DEGREES (TRANSFER APOGEE IS 111 N. MI.) 


Iteration 

*1 

— 

*2 

^3 

Gradient 3- 

Payoff 

0 

0.0 

0.0 

28.5 

0.57E-03 

0.49310948 

1 

0.167574 

0.167574 

28.1649 

0.15E-08 

0.49199596 

2 

0.172128 

0.163036 

28.1648 

0.38E-11 

0.49199583 

3 

0.171911 

0.162796 

28.1653 

0. 58E-14 

0.49199583 

4 

0.171901 

0.162804 

28.1653 

0.46E-18 

0.49199583 

5 

0.171902 

0.162805 

28.1653 

0.8 IE -26 

0.49199583 

0 

0.0 

28.5 

0.0 

0.29E-03 

0.49218163 

1 

0.0 

28.4841 

0.0158555 

0.29E-03 

0.49212881 

2 

0.169402 

28.3147 

0.0158615 

0. 91E-11 

0.49155075 

3 

0.169402 

28.3148 

0.0158285 

0.11E-17 

0.49155075 

4 

0.169402 

28.3148 

0.0158285 

0. 94E-18 

0.49155075 

5 

0.169402 

28.3148 

0.0158285 

0. HE -21 

0.49155075 

0 

28.5 

0.0 

0.0 

0.29E-03 

0.49327239 

1 

28.4835 

0.0 

0.0165287 

0.29E-03 

0.49321892 

2 

28.3165 

0.166959 

0.0165348 

0.22E-10 

0.49268099 

3 

28.3166 

0.166959 

0.0164769 

0.32E-16 

0.49268099 

4 

28.3166 

0.166959 

0.0164770 

0.27E-16 

0.49268099 

5 

28.3166 

0.166958 

— 

0.0164769 

■ 

0.21E-21 

0.49268099 
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TABLE 8. DISTRIBUTION OF PLANE CHANGES FOR IMPULSIVE 
TRANSFER BETWEEN 100 AND 110 N. MI. CIRCULAR 
ORBITS INCLINED AT 60 DEGREES (TRANSFER APOGEE IS 111 N. MI.) 


Iteration 



^3 

Gradient 3- 

Payoff 

0 

0.0 

0.0 

60.0 

0.46E-03 

1.0001054 

1 

0.0735228 

0.0735228 

59.8530 

0.50E-06 

0.99936690 

2 

0.0766450 

0. 0704655 

59. 8529 

0.26E-10 

0. 99936641 

3 

0.0765758 

0.0703931 

59.8530 

0.59E-13 

0.99936641 

4 

0.0765724 

0.0703964 

59.8530 

0.73E-16 

0.99936641 

5 

0.0765723 

0.0703963 

59. 8530 

0.81E-21 

0. 99936641 

0 

0.0 

60.0 

0.0 

0.23E-03 

0.99887366 

1 

0.0 

59.9930 

0.00698963 

0.23E-03 

0.99883858 

2 

0. 0762936 

59.9167 

0.00699231 

0.94E-10 

0.99845353 

3 

0.0762936 

59.9167 

0.006 97960 

0.34E-19 

0.99845353 

0 

60.0 

0.0 

0.0 

0.23E-03 

1.0011622 

1 

59.9929 

0.0 

0.00705916 

0.23E-03 

1.0011269 

2 

59.9221 

0.0708007 

0.00706168 

0. 15E-09 

1.0007720 

3 

59. 9222 

0.0708007 

0.00704546 

0.UE-18 

1.0007720 

4 

59.9222 

0.0708007 

0.00704546 

0. 55E-19 

1.0007720 
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TABLE 9. DISTRIBUTION OF CONSTRAINED PLANE CHANGES FOR 
IMPULSIVE TRANSFER BETWEEN 100 AND 150 N. MI. CIRCULAR 
ORBITS INCLINED AT 28.5 DEGREES (TRANSFER APOGEE IS 200 N. MI.) 


Iteration 


v 

V 

Gradient* 3 

Payoff 

Case i: 5 degree limit on 5 degree limit on 

0 

1.25 

1.25 

26.0 

0.92E-06 

0. 49218937 

1 

1,4919609 

1.2670539 

25.740985 

0.69E-08 

0. 49216474 

2 

1.4929680 

1.3367815 

25.670251 

0.23E-11 

0. 49216410 

3 

1.4934366 

1.3368224 

25.669741 

0.13E-16 

0. 49216410 

4 

1. 4934363 

1.3368255 

25.669738 

0. 12E-17 

0.49216410 

5 

i. 4934366 

1.3368256 

25. 669738 

0.40E-24 

0. 49216410 

Case 2: 1 degree limit on^, 5 degree limit on ^ 

0 

0.25 

1.25 

27.0 

0.43E-04 

0. 49432670 

i 

0. 99999878 

1.1942897 

26.305712 

0. 12E-07 

0. 49229098 

2 

0.99999969 

1.2747884 

26.225212 

0.48E-12 

0.49228996 

3 

1. 0000000 

1.2747736 

26.225226 

0. 45E-17 

0.49228996 

4 

... 

i 

1.2747719 

26.225228 

0. 12E-21 

0. 49228996 

Case 3: 5 degree limit on £ , 1 degree limit on i? 2 

0 

1.25 

0.25 

27 

0. 10E-04 

0.49291871 

1 

1.9253040 

0.78144550 

25.793250 

0. 13E-05 

0.49229192 

2 

1.4313804 

0.83753959 

26.231080 

0.57E-07 

0.49221594 

3 

1.5509740 

0.99627414 

25.952752 

0.54E-07 

0.49218585 

4 

1.4770807 

0.99641748 

26.026502 

0.27E-09 

0.49218394 

5 

1.4778344 

1.0000000 

26.022166 

0.55E-12 

0.49218345 

6 

1.4776119 

1 

26.022388 

0.11E-18 

0.49218345 

7 

1.4776123 

1 

26.022388 

0.15E-17 

0.49218345 

8 

1.4776120 

1 

26.022388 

0.35E-24 

0.49218345 

Case 4: 1 degree limit on^, 1 degree limit on 

0 

0.25 

0.25 

28.0 

0.52E-04 

0. 49502071 

1 

0. 98877934 

0.76936973 

26.741851 

0.97E-07 

0.49236042 

2 

0.99685410 

0. 98658947 

26. 516556 

0. 57E-08 

0. 49230787 

3 

0. 99990077 

0. 99604523 

26.504054 

0.38E-09 

0. 49230480 

4 

0.99991152 

0.99849260 

26.501596 

0.22E-09 

0.49230450 

5 

0. 99999996 

0.99999993 

26.500000 

0.71E-13 

0.49230427 

6 

1.0000000 

1.0000000 

26.500000 

0.19E-18 

0. 49230427 

7 

1.0000000 

1.0000000 

26.50000 

0.61E-19 

0.49230427 

8 

1 

1 

26.5 

0.21E-27 

0.49230427 


a. Initially, x^ = Xg = 30 degrees 


b. Gradient = 
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